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Electron systems that possess lig:ht-like dispersion relations or the conical 
Dirac spectrum, such as graphene nj, |2n and bismuth , have recently been 

shown to harbor unusual collective states in high magnetic fields. Such states 
are possible because their light-like electrons come in spin pairs that are chi- 
ral, which means that their direction of propagation is tied to a quantity called 
pseudospin that describes their location in the crystal lattice. An emerging di- 
rection in quantum materials research is the manipulation of atomic spin-orbit 
coupling to simulate the effect of a spin dependent magnetic field, in attempt 
to realize novel spin phases of matter This effect has been proposed to 

realize systems consisting of unpaired Dirac cones that are helical, meaning their 
direction of propagation is tied to the electron spin itself, which are forbidden to 
exist in ordinary Dirac materials such as graphene or bismuth. Such unpaired 
spin-helical Dirac cones are predicted to exist at the edges of certain types of 
topologically ordered materials M, yi and support exotic collective states that 

n n n n 

exhibit non-Maxwell electromagnetic behavior (91. IIQ. II IL Il2[]. charge fractional- 
ization |13|] and topological proximity effects when interfaced with conventional 
electron systems |14l. Il5[|. The experimental existence of topological order cannot 
be determined without spin-resolved measurements {3, Q, [l6|- Here we report 
a spin- and angle-resolved photoemission study of the hexagonal surface of the 
Bi2_a:Mna^Te3 series, which is found to exhibit a single helical Dirac cone that is 
fully spin-polarized. Our observations of a gap in the bulk spin-degenerate band 
and a spin-resolved surface Dirac node close to the chemical potential show that 
the low energy dynamics of Bi2Te3 is dominated by the unpaired spin-helical 
Dirac modes that are entangled over the macroscopic dimensions of the crystal. 
Our spin-texture measurements prove the existence of a rare topological phase 
in this materials class for the first time, and makes it suitable for novel 2D Dirac 
spin device applications beyond the chiral variety or traditional graphene. 

Three-dimensional topologically ordered materials represent a new phase of quantum 
matter where a combination of relativistic and quantum entanglement effects dominate 51. 
They were recently proposed p, 116|] and shortly afterwards first discovered in the Bii_a;Sba; 
alloy material [17]. In these systems, spin-orbit coupling gives rise to non-current carrying 
electronic states in the bulk and robust conducting states along the edges. In contrast to 
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graphene, which has four Dirac cones (2 doubly degenerate cones at the K and points 
in momentum space) [2], the remarkable property of topological edge states is that their 
dispersion is characterized by an odd number of spin-polarized Dirac cones. Therefore, it 
would be impossible to prove the existence of topological order in materials without probing 
the spin degrees of freedom as demonstrated recently by Hsieh et.al.l^]. Odd spin edge 
metals are expected to exhibit a host of unconventional properties including a 
(half-integer) quantum Hall effect [t^ and immunity to Anderson localization {3, 2lj]. 
most exotic physics, however, may occur at the interface between a topologically ordered 
material and an ordinary ferromagnet or superconductor, where electromagnetic responses 
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i;hat defy Maxwell's equations [9|, llO, Uj| and excitations that obey non-Abelian statistics 
1^ . [isl have been predicted. 

The key to testing such experimental proposals is to find the most elementary form of 
helical Dirac spin states that are fully spin polarized. Despite an active search no spin 



helical Dirac cones have been experimentally observed in Bii_a;Sba; |17|, |25| . Bi9Se3 jl8l| or 
Bi2Te3[l9|. At present no spin-resolved surface band calculation on Bi2Te3 exists and spin- 
texture is theoretically unknown. Here we report both a conventional and spin-resolved 
ARPES investigation of single crystals of Bi2Te3 and Bi2_xMna;Te3. Remarkably, a single 
helical Dirac cone is realized on the (111) surface, with fully polarized spins rotating by 
360° around a single Fermi surface. Although we observe a time dependent band bending 
of the surface states, a long standing problem with Bi2Te3 [19], we show a method that 
the characteristic relaxation rate can be systematically controlled by adding small amounts 
of Mn, and the spin Dirac node can be fairly stabilized near the Fermi level E^. Using 
a synchrotron light source with a tunable photon energy (hu)^ we show that the bulk-like 
states of Bi2Te3 are insulating with the valence band maximum lying around 100 meV below 
E'i?, consistent with the long-theorized six-peak dispersion model. Fully gapped bulk-like 
states together with a spin-polarized Dirac node lying at a very low energy makes Bi2Te3 
and slightly Mn-doped Bi2_xMna^Te3 suitable for providing 2D helical spin-textured Dirac 
fermions for spin device integration. 

High resolution ARPES data were taken at Beamlines 12.0.1 and 7.0.1 of the Advanced 
Light Source in Lawrence Berkeley National Laboratory. Spin-resolved ARPES data were 
taken with the COPHEE spectrometer [22^ at the SIS beamline of the Swiss Light Source 
using horizontally polarized 22 eV photons. At a pass energy of 3 eV, the energy resolution 
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for spin-detection was set at around 30 meV. The COPHEE spectrometer is capable of 
measuring all three spatial components of the spin polarization vector for any point in 
reciprocal space from which a spin-resolved band structure is directly imaged. Data were 
taken from single crystal Bi2_a:Mna;Te3 cleaved along its (111) surface (Fig.[l](a)) in ultra high 
vacuum (UHV) at pressures better than 5 x 10~^^ torr and maintained at a temperature 
of 15 K. Single crystals of Bi2_xMna^Te3 were grown by melting stoichiometric mixtures 
of elemental Bi (99.999 %), Te (99.999 %) and Mn (99.95 %) at 800°C overnight in a 
sealed vacuum quartz tube. The crystalline sample was cooled over a period of two days 
to 550°C, and maintained at the temperature for 5 days. It was then furnace cooled to 
room temperature. The crystal structure of Bi2_a:Mna;Te3 is rhombohedral with space group 
Dl^{R3m). It can be visualized as a stack of hexagonal atomic layers, each consisting of only 
Bi/Mn or Te. Five atomic layers are stacked in a close-packed fee fashion along the [111] 
direction in order Te(l)-Bi-Te(2)-Bi-Te(l), in a quintuple layer, and cleavage takes place 
naturally between such layers. The theoretical calculations were performed with the LAPW 



method in slab geometry using the WIEN2K package [23|]. GGA of Perdew, Burke, and 
Ernzerhof [24] was used to describe the exchange-correlation potential. Spin-orbit coupling 
was included as a second variational step using scalar-relativistic eigenfunctions as basis. 
The surface was simulated by placing a slab of six quintuple layers in vacuum. A grid of 
35x35x1 points were used in the calculations, equivalent to 120 k-points in the irreducible 
BZ and 2450 k-points in the first BZ. 

The simplest signature of topological Z2 order is the existence of a single spin-polarized 
Dirac cone on its surface, with the Dirac node located at a momentum k^^ in the surface 
Brillouin zone (BZ), where satisfies k^ = -k^ + G, G being a surface reciprocal lattice 
vector [7]. Because this node is protected by time-reversal-symmetry 3], the spin polariza- 
tion is guaranteed to rotate by 360° about its Fermi surface [25|] (Fig. [I](b)). Figure [I](c) 
shows a spin-integrated ARPES intensity map of Bi2Te3 (111) taken at a constant energy 
which spans several surface BZs. The density of states at is clearly localized only 
to a small area around the k^ = F points, suggesting a very small concentration of surface 
carriers. A high resolution intensity map around F shows that the density of states is dis- 
tributed about a ring enclosing F (Fig. [I](d)). This allows for a calculation of the surface 
carrier density n using the Luttinger theorem n = 'Kk']^/ {Abz^r)) where kp the Fermi 
momentum, Abz is the area of the surface BZ and Aji is the area of the surface unit cell. 
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which yields n ^ 1.7 x 10^^ cm~^. The dispersion of the Dirac cone (Figure [T](e)), imaged 
an hour after cleavage, shows a Dirac node located at a binding energy (Eb) of around -130 
meV. Because experimental proposals to measure non-Maxwell electrodynamics H, 26] 



or non-Abelian quasiparticles Q require the surface state to be made insulating by 
lifting the spin degeneracy of the Dirac node with a magnetic field, a Dirac node close to 
means that this can be achieved using smaller magnetic fields. 

ARPES spectra taken as a function of time show that the binding energy of the Bi2Te3 
surface Dirac node exhibits a pronounced time dependence, increasing from Eb ^ -100 meV 
8 minutes after cleavage to Eb ^ -130 meV at 40 minutes (Fig. [2](d)-(f)), in agreement with 
a previous report Such behavior has been attributed to a downward band bending 

near the surface (Fig. [21(b)) that is caused by the breaking of inter-quintuple layer van der 
Waals Te(l)-Te(l) bonds (Fig. [2](a)) and the formation of Te(l)-Bi bonds [271, bsl, which 
increases the electron density near the surface. To gain insight into the surface state forma- 
tion process, we carried out ab initio calculations with optimized lattice parameters. The 
calculations (Fig. [2](c)) show that the surface Dirac cone only appears when spin-orbit cou- 
pling is switched on and that its dispersion closely matches the late-time ARPES spectrum 
(Fig. [2](f)). The calculated position of the Dirac node lies in the bulk band gap, which 
corroborates our experimental data showing a strong intensity near the Dirac node and a 
drastic weakening away from F as the surface band merges with the bulk bands to form 



resonance states 
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2Mi. The slow dynamics of the band bending process suggests that 



the local surface lattice structure is coupled to the local surface charge density. Equilibration 
involving both atomic and charge redistribution may thus be significantly delayed by struc- 
tural defects such as site inversions or vacancies. By systematically increasing the defect 
concentration through Mn/Bi substitution, we show that band bending can be slowed by up 
to 10 fold (Fig. [2](g)-(i)), allowing a wider range of its intrinsic relaxation time scale to be 
accessed. ARPES valence band spectra (Fig. [2](k)) of Bii.95Mno.o5Te3 taken over a 15 hour 
period show that the positions of the valence band edges shift downward by a total energy 
of around 100 meV (Fig. [2](k)), representing the total magnitude of band bending A. 

Our theoretical calculations show that the bulk valence band maxima (VBM) in Bi2Te3 
lie in the TZL planes of the bulk BZ around 15 meV below E^ (Fig. [3](b)), giving rise to 
VBM in each of the six mirror planes in agreement with theoretical proposals 
order to measure the locations of these bulk VBM in the three-dimensional BZ, we performed 
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a series of ARPES scans along the trajectories shown by the red hnes in Figure [3](a) with 
varying incident photon energies to access multiple values of fc^. In order to identify subtle 
changes in the band structure with changing /c^, it is desirable to work with a band structure 
that is stable in time. Therefore all /I'U-dependent scans were taken more than an hour after 
cleavage. However, because ARPES is only sensitive to the topmost quintuple layer at 
the photon energies we use [29], the three-dimensional band structure we measure is only 
bulk-like in the sense that the energy positions may differ from the true bulk states due 
to near-surface band bending effects. Figures [3](c)-(e) show a series of ARPES intensity 
images of the band dispersions along momentum space trajectories in the kx-kz plane taken 
using photon energies of 31 eV, 35 eV and 38 eV respectively (Fig. [3]). The Dirac cone near 
Ejp centered at kx = shows no dispersion with photon energy, supporting its surface 
state origin. In contrast, a strongly kz dispersive hole-like band is observed near kx = 0.27 
A~-^, whose peak is found to lie closest to at hu = 35 eV (Fig. [3](d)). Using the free 



electron final state approximation 
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29|, this VBM lies close to the (0.27, 0, 0.27) point 



in the bulk BZ, which agrees well with our calculations. The energy position S of the VBM 
lies approximately 100 meV below based on its peak position in the energy distribution 
curves (Fig. [3](g)), which is lower than that found in our calculations by ^^85 meV. This 
discrepancy agrees well with magnitude of the near-surface band bending A that is also 
responsible for the lowering of the surface band energies (Fig. [2j). We conclude that the 
three-dimensional states near the surface region behave as a large band gap insulator. 

To determine the helicity of the surface Dirac cones, we performed spin-resolved ARPES 
scans jsil using a double Mott detector configuration (Fig. [4](d)), which allows all three 
components of the spin vector of a photoelectron to be measured [32]. We analyzed the 
spin-polarization of photoelectrons emitted at Eb = -20 meV along momentum space cuts 
traversing the surface states (Fig. H]). Although scans were taken hours after cleavage in 
vacuum to allow time for the surface states to stabilize, the surface band dispersion continued 
to slowly evolve over the course of the scans. In effort to compromise between the effects 
of changes in the surface state dispersion with time, increasing surface contamination over 
time, and the typical hours-long counting times needed for spin-resolved measurements, we 
collected data with a level of statistics sufficient to show the spin-polarized character of the 
surface states. Figures [11(a) and (b) show the spin polarization spectra of the x, y and z 
(out-of-plane) components along the F-M direction. In the x and z directions, no clear signal 
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can be discerned within the margins of statistical error. In the y direction on the other hand, 
clear polarization signals of equal magnitude and opposite sign are observed for the surface 
bands on either side of F, evidence that they form a time-reversal-invariant spin-orbit split 
pair. Because F-M is along a line of mirror symmetry in the Bi2Te3 crystal, it is expected 
that the spins are constrained along the ±ii direction along this cut. A fit based on the 
two-step routine developed by Meier et al. [32] is consistent with 100% polarized spins that 
are directed completely along the ±y direction as shown in Figure [4](e) . Figure [4](c) shows 
the spin-resolved spectra (/J'^) calculated from Py according to /J = Itot{^ + Py)/"^ ^^id 
I^y = Itot{l - Py)/2, where Ifot is the spin-averaged intensity, which shows a peak to peak 
separation in agreement with 2kF. The non-Lorentzian lineshape of the /J and curves 
and their non-exact merger at large \kx\ is likely due to the time evolution of the surface 
band dispersion. 

Our combination of conventional and spin-resolved ARPES results show that 
Bi2-xMna^Te3 hosts a single unpaired helical Dirac cone on its (111) surface that is fully 
spin-polarized, opening rare opportunities for Dirac spin materials research beyond the 
chiral variety for the first time. Our finding of large-gap insulating behavior of the near 
surface bulk-like spin-degenerate states in the x = compound, together with a low lying 
spin-resolved Dirac node makes Bi2Te3 suitable for studying topological effects via interface 
transport. 
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kx (A^) Eb (eV) 

FIG. 1: . A single massless Dirac cone on the surface of Bi2Te3. (a) Photograph of the 
cleaved (111) surface of Bi2Te3. (b) Schematic of the (111) surface Brihouin zone with the four 
time-reversal-invariant momenta (r,3xM) marked by blue circles. A single spin-polarized Fermi 
surface enclosing F that arises from a spin-textured Dirac cone is the direct signature of topological 
order, (c) Constant energy ARPES intensity map of Bi2Te3 collected at Ejr using hv = 100 eV. 
(d) High resolution view of the Fermi surface, using hu = 35 eV, over a momentum space window 
demarcated by the white box in (c). (e) ARPES intensity map of the gapless surface state bands. 
The red dotted lines are guides to the eye showing the Dirac dispersion, (f) Energy distribution 
curves of the data shown in (e). 
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FIG. 2: . Slow dynamics of the surface Dirac cone dispersion, (a) The crystal structure of 
Bi2Te3 viewed parallel to the quintuple layers. The Te(l) 5p orbitals that form the inter-quintuple 
layer van der Waals bonds are shown in yellow, (b) Schematic of the band bending of the bulk 
valence band maximum (VBM) near the cleaved surface, (c) Calculated band structure along the 
K-f-M cut of the Bi2Te3(lll) BZ. Bulk band projections are represented by the shaded areas. The 
band structure results with spin-orbit coupling (SOC) are presented in purple and that without 
SOC in yellow. No pure surface band is observed within the bulk band gap without SOC (black 
lines). One pure gapless surface band crossing is observed when SOC is included (red lines). 
Inset shows enlargement of low energy region near F. Spin-integrated ARPES spectra of Bi2Te3 
along the f-M direction taken with hu = 30 eV (d) 8 mins, (e) 20 mins, and (f) 40 mins after 
cleavage in UHV. Analogous ARPES spectra for Bii gsMno.osTca (g) 15 mins, (h) 4 hours and (i) 
9 hours after cleavage, showing a slower relaxation rate. Red lines are guides to the eye. (j) High 
energy valence band ARPES scan of Bii.95Mno.o5Te3 and (k) the energy distribution curves at F 
taken at different times. 
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FIG. 3: . Observation of insulating bulk-like states supporting a six-peak electronic 
structure, (a) Bulk rhomb ohedral Brillouin zone of Bi2Te3. According to LDA band structure 
calculations [sj], six valence band maxima are located at the b points that are related to one another 
by 60° rotations about z. The red hues show the momentum space trajectories of the ARPES scans 
taken using hv — 31 eV, hv — 35 eV and hv — 38 eV. (b) Calculated band structure through the b 
point along cut 2 that has been rigidly shifted to match the superimposed corresponding ARPES 
cut to account for the Fermi level shift in real materials, (c) to (e) show ARPES intensity maps 
along the hv — 31 eV, 35 eV and 38 eV trajectories respectively, obtained one hour after sample 
cleavage. The in-plane momentum component of the b point is marked by an arrow in (d), and 
the energy of the valence band maximum relative to Ep is denoted 5. (f) to (g) show the energy 
distribution curves corresponding to images (c) to (e) respectively. 
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FIG. 4: . Helicity of surface Dirac cone with a non-trivial topology, (a) Measured y 
component of spin-polarization along the F-M direction at Eb = -20 meV, which only cuts through 
the surface states. Inset shows energy position of this cut together with fitted spin-polarization 
directions. The sohd hne is a fit based on the procedure of Meier et al js^]. (b) Measured x 
(orange diamonds) and z (green circles) components of spin-polarization along the f-M direction 
at Eb = -20 meV. (c) Spin-resolved spectra obtained from y component spin polarization data, (d) 
Geometry of the Mott detector axes, (e) Fitted values of the spin polarization vector demonstrating 
helicity of the Dirac cone. The largest source of uncertainty in the fitted angles comes from the 
time dependence of the surface state spectrum, (f) Spin-integrated ARPES intensity map of the 
surface state Fermi surface with its the measured helicity superimposed. This originates from the 
spin-polarized Dirac cone shown to its right. 
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